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Abstract 

Top-quark physics plays an important role at hadron colliders such as 
the Tevatron at Fermilab or the LHC at CERN. Given the planned 
precision at these colliders, precise theoretical predictions are re- 
quired. In this paper we present the complete electroweak correc- 
tions to QCD-induced top-quark pair production in quark-antiquark 
annihilation. In particular we provide compact analytic expressions 
for the differential partonic cross section, which will be useful for 
further theoretical investigations. 



I. Introduction 



At ongoing and upcoming collider experiments, top-quark physics will play a cen- 
tral role. Although the top-quark was discovered already 10 years ago, direct mea- 
surements of its properties are still rather limited. In particular most of the quantum 
numbers are only constrained from indirect measurements such as the electroweak 
precision observables. In the near future the hadron colliders Tevatron at Fermilab and 
LHC at CERN will provide unique possibilities for detailed measurements in the top 
sector. A necessary requirement for these analyses is the precise theoretical under- 
standing of reactions involving top-quarks. At hadron colliders both single top-quark 
production as well as top-quark pair production have been studied extensively in the 
past. The differential cross section for top-quark pair production is known to next-to- 
leading order (NLO) accuracy in quantum chromodynamics (QCD) ||TJ|5J|51I11I21|. In 
addition, the resummation of logarithmic enhanced contributions has been studied in 
detail in Refs. [71 [HI El HOI El- Recently also the spin correlations between top-quark 
and antitop-quark were calculated at NLO in QCD [E1[T2|. In Refs. 021 [ll] the elec- 
troweak corrections were investigated. However, for the quark-antiquark annihilation 
process, only the electroweak vertex corrections were considered — the contributions 
from box diagrams were ignored. It is well known that in the high energy region 
s ^> mw.z the weak corrections can be enhanced by the presence of large logarithms 
(see e.g. Refs. lfl"5l IT6ll and references therein) which justifies a detailed study of all 
contributions. More recently the electroweak corrections for £-quark production were 
re-analysed in Refs. ifTTlfTFIl . In particular, it was again confirmed that the weak cor- 
rections can lead to sizable corrections for specific observables. For a more detailed 
theoretical investigation of these effects, it is useful to have short analytic expressions 
available. The aim of this work is to recalculate the weak corrections to top-quark pair 
production — including the contribtution from box diagrams — and give compact an- 
alytic results. Note that we do not consider here the pure photonic corrections, which 
form a separate gauge-invariant subset. 

The outline of the paper is as follows. In section 2 we present the calculation of 
the virtual electroweak corrections to top-quark pair production in quark-antiquark 
annihilation. The contributions involving box diagrams are infrared-divergent. The 
singularities cancel when the virtual corrections are combined with the corresponding 
real corrections, which we calculate in section 3. In section 4 we discuss some checks 
we performed and give numerical results for the total cross section. 



II. Virtual corrections 

In this section we present the calculation of the electroweak corrections. We work in 
the 't Hooft gauge (i?^-gauge) with the gauge parameters ^' set to 1. In this gauge, 
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apart from the physical fields, also unphysical fields contribute. In particular we have 
to consider the contribution from the fields denoted by %, which are related to the 
longitudinal degrees of freedom of the gauge bosons. In principle in the i?^-gauge also 
ghosts need to be considered to cancel unphysical degrees of freedom. To the order 
where we are working, the ghosts do not contribute. In addition, given that we neglect 
the masses of the u,d,c,s quarks the unphysical fields (|) and % only contribute in the 
vertex corrections to the final gluon-top-antitop vertex. The renormalization is done in 
renormalized perturbation theory. That is the bare Lagrangian L is rewritten in terms 
of renormalized fields and couplings: 



^(^0,^0,^0,^0) 



= L (Z^ 2x ¥ R , z\ /2 A R , Z m m R , Z g g R ) 

= L( x V R ,A R ,mR 1 g R ) + L ct ( x i , R,AR,mR,g R ). (II. 1) 



The contribution L{^ R ,A R ,m R ,g R ) gives just the ordinary Feynman rules, but with 
the bare couplings replaced by the renormalized ones. Some sample diagrams are 
shown in Fig. III. 11 The complete list of Feynman rules can be found for example in 
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Figure II. 1: Sample diagrams for the virtual corrections. 



Ref. |[T91l. The second contribution in Eq. ( IH.1I) L ct (*¥ R ,A R , m Rl g R ) yields the counter- 
terms, which render the calculation ultraviolet (UV)-finite. The diagrams needed here 
are shown in Fig. 111.21 Note that although the electroweak corrections appear here 
in one-loop approximation, they are the leading-order electroweak contribution. The 
interference term of the amplitude M (qq — > y, Z — > tt) with the corresponding QCD 
amplitude vanishes as a consequence of the specific colour structure. Terms of order 
GCyCC are therefore absent. Thus no renormalization of the coupling constants is required 
at the order under consideration here. This is different from an electroweak correction 
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Figure II. 2: Counterterm diagrams. 



to an electroweak amplitude, which would not be UV-finite without coupling-constant 
renormalization. The whole contribution from the renormalization is given by: 

h\M | 2 = 2(dZ q + dZ t )\M qq ^ ti \ 2 , (11.2) 

where Z q , Z t denote the wave-function renormalization constants of the incoming light 
quark and the outgoing top-quark (Z; = 1 + 8Z ; ). The squared leading-order QCD 
amplitude \M qq -^ t f\ in d dimensions is given by: 

l^^fl 2 = l6n 2 a 2 {N 2 - 1)(2-(3 2 (1 -z 2 ) -2e), (II.3) 

where N is the number of colours, a s the strong coupling constant and [3 the velocity 
of the top-quark in the partonic centre-of-mass system: 



13 = ^1-4^ (H.4) 

(s denotes the partonic centre-of-mass energy squared). The cosine of the scattering 
angle is denoted by z- The parameter of dimensional regularization e is defined by 

d = 4 - 2e. (H.5) 

For the renormalization of the quark fields we use the on-shell scheme. The renor- 
malization constants in this scheme in terms of self-energy integrals and derivatives 
thereof can be found for example in Ref. [ 19 1. Before presenting the results, let us add 
a few technical remarks. We used the Passarino-Veltman reduction scheme [|2H to 
reduce the tensor integrals to scalar one-loop integrals. For the scalar integrals we use 
the following convention: 



in 2 J (£ 2 - mf + ie) • • • 



= ^2 d d £ - 2 \J J - ■ (H.6) 



For the UV-divergent integrals we define the finite part for the one-point integrals Ao 
and the two-point integrals Bq through 

A (m 2 ) = m 2 A+A (m 2 ), 
B (p 2 ,mlm 2 2 ) = A + B (p 2 ,m?,m|), (H.7) 
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with 

A = (47i) e r(l + e)i = ^ -y+ln(47c) + 0(e). (II.8) 

The vertex corrections do not contain infrared or mass singularities (IR singularities). 
They contain only UV singularities, which are removed by the aforementioned renor- 
malization. On the other hand the contribution involving box diagrams are UV-finite 
but contain IR singularities. In order to regularize the IR singularities, we use dimen- 
sional regularization. To simplify their determination, we express the J-dimensional 
four-point scalar integrals Dq in terms of the (d + 2) -dimensional four-point integrals 
Dq +2 and a combination of three-point integrals in d dimensions. This can be done by 
the following relation 

D d Q +z = -2%D d 21 {11.9) 
where the box integral in 6 dimensions is defined by 

D '~ {PuP2,P3>Pi-P2>P2-P3,Pl-P3,mi ,m 2 ,m 3 ,m 4 )[l,^,£/ v , . . .] = 
1 f d 6 £ ( [^Wv,--\ 



in 2 J I (e 2 -mi 2 + ie)((e + p l ) 2 -m 2 2 + ie) 

X ((£ + pi+ p 2 ) 2 - m 3 2 + ie){{£ + P1+P2+ P2,) 2 - m A 2 + ie) } ^' 10 ^ 

and £>2 7 is the coefficient of the metric tensor appearing in the Passarino-Veltman 
decomposition EUll of the tensor integral 

D d =%£ v ] (11.11) 

(see Eq. (F.3) in Ref. Il20l0 . which in turn can be expressed as a linear combination of 
the scalar box integral Dq and scalar triangle integrals Cq in d dimensions. Owing to 
the finiteness of the box integrals in 6 dimensions the infrared singularities appear only 
in the three-point integrals. 

For the presentation of the results it is convenient to use the leading-order QCD cross 
section which is given by: 

^— =a (2-(3 2 + |3¥) (11.12) 

dz 

with 

1 AJ 2 — i R 

° 0= 8 m '^V^? (IU3) 
A factor 1 / {AN 2 ) from averaging over the incoming spins and colour is included. In 
a first step we present the corrections to the cross section from the contribution of Z 
boson and W boson exchange to the light quark-gluon-vertex. The form of this vertex 
remains unchanged, its normalization is shifted by a factor (1 +8F^ n W+Z ), which 
leads to a shift by 

j.-rlrL.W+Z Born 

^— = ^2Re5^ n < w+z (11.14) 

dz dz 
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with 



6F, 



In,W+Z 
D 



1 a 

871 



(gf +gf)MPz) + 2gw 2 Mp w ) 



and 



fi(x) = 1+2 (l+ln(x))(2z+3)-2(l+x) z Li 2 1 + - - 



where we used the definition 



Pi 



m, 



(11.15) 



(11.16) 



(11.17) 



^ 2 

In Eq. (III. 151) . a = |^ denotes the electromagnetic coupling. 

The Cabibbo -Kobayashi-Maskawa mixing matrix has been set to 1. The vector and 
axial vector couplings of neutral and charged currents are given by 
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gw = 



2s\yCw 
1 

2s\yC iy 
1 

2y/2sw 



[T«-2s w 2 Q q ), 



J 3 ' 



(11.18) 
(11.19) 
(11.20) 



where Q q describes the electric charge in units of the elementary charge e, sw is the 
sine of the Weinberg angle (s\y = shifty) , cw = cos(tV)X and T q denotes the weak 
isospin. For the vertex corrections to the final vertex we split the result into the contri- 
bution from Z boson exchange, W boson exchange, Higgs exchange and the contribu- 
tions from the would-be Goldstone bosons including the respective counter terms: 



do Vm - = do Fin ' z + do ¥m - w + do Vm ' H + do Fin * + do Fm -*. 



(11.21) 



Note that only the sum has a physically meaningful interpretation. For the individual 
contributions we obtain 



do 



Fin.,Z 



a f 



o ^j-2(^ 2 +^ 2 )(l+^) 
(g[ 2 + g t a 2 )(M'nz 2 )-M>nt 2 )) 



+ 8 
- 4 



dz 

l+z 2 

(1-P 2 > 
1 



(1-P 2 )P 2 
+ fl-5z 2 -3(3 2 (l-z 2 



(l - (3 2 + 2f3 4 ( 1 - z 2 ) + 7(3 V - 3z 2 ) p z {g t2 + g' 2 ) 



^ 2 -(3 + z 2 -p 2 (l-r 



flj(3,4) 



+ 2 



2\„f2 



p2/ 2 (z,(3)p z (^ 2 + ^ 2 ) -/2(z,P)(^ 2 -3^ 2 ) -4(1 +z 2 K 



4(1,3) 
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+ 2s 



4(/ 2 feP)^ 2 - (l +z 2 )g t v 2 )p z + ^f 2 (z^)p 2 z (g t v 2 + g t a 2 ) 



(l + p 2 )(2-p 2 (l-z 2 )V v 2 -(2-3p 2 + 5(3V + 3p 4 (l-z 2 ))^ 



r 3 



+ 2s 



2p z (^ 2 +^ 2 ) + (l-p 2 )(^ 2 -3^ 2 ) 



with 



/2 (z,P) = l-3z 2 -2p 2 (l-z 2 ). 



(11.22) 



(11.23) 



The integrals are defined in the appendix. For the contribution from the W boson we 
obtain: 



do 



Fin.W 



a 



dz 



°o^gw z i -2(l+z / ) + 8 



l+z 2 



(1 - P 2 )* 



(A (m w 2 )-A (m b 2 )) 



(l - 3z 2 - p 2 + 7pV + 2p 4 (l - z 2 )) (p. - p b ) 



p 2 L(l-p 2 ) 
+l-3z 2 -5p 2 (l+z 2 )-2p 4 (l-z 2 )]fl 4 ) (l,2) 



1 

+ pi 

s 

+ 



4/ 2 (z,p)(p w -p,) + l-3z 2 -5p 2 (l+z 2 )-2p 4 (l-z 2 ) 55(1,3) 



4P 2 L 



8/ 2 (z, p) (2(p w - p b f - (1 + p 2 )p 6 ) + 1 - 3z z - 16p 2 - 4P 2 z 2 + P 4 



llpV-2p 6 (l-z 2 )+8(l-3z 2 -2p 4 (l-z 2 )-3p 2 (l+z 2 ))p w ]c^ 



+ ,(2-p 2 (l-z 2 )) 4(p w -p,)-(l-p 2 ) ^^(1,2)| , 2 



\p*=m t * 



(11.24) 



The contribution from the Higgs boson is given by 

-0+; 2 ) 



do Fin.,H a m 2 



+ 4 



l+z 2 - - 

(A (ra// ) -A (m f 2 )) 



(1-P 2 > 



^ ( 1 - p 2 + 2P 4 ( 1 - z 2 ) + 7P 2 z 2 - 3z 2 ) Ph 



(l-p2)P 
+ 2(l-z 2 )(l-p 2 )l^(l,2) 





1 - 


+ 


pi 




2s - 


+ 


pi 



2/ 2 (z, p)p„ + p 2 (5 - 3z 2 - 4p 2 (l - z 2 ))] Sq(1 , 3) 

/ 2 (z,p)p 2 ; + 3p 2 (l-z 2 )(l-p 2 )p w -p 2 (l-p 2 )(2-p 2 (l-z 2 ))]c ( 



+ 2s (2 - P 2 ( 1 - z 2 ) ) [p H - ( 1 - P 2 )] B l ( 1 . 2 ) 



p 2 =m t 2 



(11.25) 



For the unphysical fields % and (]) we find 



da Fin -'X 



2n Sa mz 2 



2 r 1.-2 

2(1 -: 2 ) 8- 1 ~ - 



a "' 2m/ J ^ ' (AoK 2 )-Ao(m z 2 )) 



(1-pV 



+ 



P 2 



2/ 2 (z,(3)p z + (3 2 (l+z 2 )l5^(l,3) 



1 - p z + 2p 4 (l - z l ) + 7(3 V - 3z 2 pX(3, 4) 



45 
+ P2 



p 2 (l-p 2 ) 

/ 2 (z,p)p 2 + p 2 (l-z 2 )(l-p 2 )p 



+ 4,(2-(3 2 (l-z 2 ))p z -^^(3,4)| p2 _ (2 }, 



(11.26) 



- l 6 (l-3z 2 -p 2 + 7p 2 z 2 + 2p 4 (l-z 2 ))p^-16p 2 (l-p 2 ) 2 (l-z 2 )^ 

+ (1 -p 2 ) 2 (l -3z 2 + 3p 2 (l +z 2 ) -2P 4 (1 -z 2 ))p; 1 ]fi5(l,2) 
1 



+ 



8p 2 



4/ 2 (z,(3)(j.v-4^p,-2(3 2 ^) 



+ (1 -p 2 )(l -3z 2 + 3p 2 (l +z 2 ) -2(3 4 (1 -z 2 )) Pw 1 l^(l,3) 



32p 2 



16/ 2 (z,(3)(j.vPw + 4p 2 ^-8p 2 )+8(l-(3 2 ) 2 (l-3z 2 + 2(3 2 (l-z 2 )) 

- 16 (l - 3z 2 + p 2 + 9pV + 2p 4 (l - z 2 )) p^ 

- 4( 1 - p 2 ) (l - 3z 2 - 1 1 p 2 - 3p 2 z 2 + 2p 4 ( 1 - z 2 ))^ 



+ (1-p 2 ) 2 l-3z 2 -7p 2 + p 2 z 2 + 2p 4 (l-z 



r 4 



+ 8 



4j,-16p^ + 8(l-p 2 )^- 



(1-P 



2\2. 



x (2 - - z 2 )) ^^(1,2) | p2 ^ 2 } , 



(11.27) 



where we used the abbreviations: 



y s = l-p 2 + 4p fc , 



(11.28) 



m b 2 

)>= 2- 

mw 



(11.29) 



Again only terms proportional to g{ 2 ,g t2 or gw 2 are present. Let us now discuss 
the contribution from the box diagrams. To the order aOj 2 considered here, we can 
distinguish two different contributions: 

1 . The (box-type) electroweak correction to the QCD Born amplitude, interfering 
with the QCD Born amplitude. 

2. The QCD box diagram interfering with the electroweak Born amplitude. 



In the following we will call the first the electroweak-box (EW-box), and the second 
the QCD-box contribution. For the EW-box we obtain 



a 

= G - 



+ 



Ja EW-box 

2(l-2p 2 + p 2 z 2 ) 



2(1 -p 2 ) 



7t\ p(l-p 2 z 2 ) 



z(l + 2p 2 -pV)^ v + 2p^ 



*i(M) 



l-Pz 
2(1~P 2 ) 

P(i-pV) 

2 -(l-2p 2 + pV 



^1(2,4) + ^ 



z(l+2p 2 -pVW v + 2^/p 



flo(l,3) 



*1(3,4) 



1+pz 



5V6v SaSa 



+ s- 



2(1 -P 2 



5o(2,4) 
; (i i-(3V-p¥-p 4 + 2p 2 



P(i-pV)(i- Pz 
+ (pV + 2pV - 2 - p 4 ) Pz - (-2P 2 - 1 + p¥) P 2 )^ v 



+ 

+ s 



2p(pV-(i-pV) Pz + p 2 )^ 



C (l,3,4) 



(2(pV-2p 2 + l)p 



(1-P^)(l-Pz) 2 

- (4 - 5p 2 + 2p 3 z + p 2 z 2 - 2p V - p 4 + p V) p 2 
+ p(p 3 z 4 + 4z- 2p V - 4p + p 3 ) p z + 3pV - 2p V 

+ 2-4pz-p 2 -pV + pV + 2p 3 z)^t 
+ 2 ((PV - 2p 2 + l)p 3 + (-3 + 4p 2 - 2p 3 z + pV - pV + p z )p 
+ (3 - 3pz + p 3 z - 3p 2 + 2p 2 z 2 )p z - 2pV + p 2 - p 3 z + p V 
+ 2Pz-l)«&' fl ]cj(2,3,4) 
1 



— s 



(1 + P^)(l-Pz) 2 



(3P 2 z 2 + 2p V + 2 + 4pz - p 2 - p V 



+ pV - 2p j z + p(p 3 z 4 - 4z - 2(3 V - 4p + p j )p z 
- (4 - 5(3 2 - 2(3 3 z + p 2 z 2 + 2|3 V - (3 4 + (3 V) p 2 

+ 2(p 2 z 2 - 2(3 2 + 1 ) p 3 ) gtgl + 2 ( + 2(3V - p 2 - p 3 z + (3 V + 2pz + 1 



+ (-3-3pz + p 3 z + 3p 2 -2p 2 z 2 )p 



+ (3-4p 2 -2(3 3 z + (3V + (3V + (3z)p 



{2 2 



- (pV-2p 2 + l)p|W fl 



C 2 (2,3,4) 



+ 



2P 2 (1 



(1-Pz)(l-Pz 



3 - p 2 - 2M 1 - Pz) - ( 1 - 2|3z + P 2 ) p z + 2p 2 gU 



+ 2(p z -(l-pz)p z + p 2 )^ 
2p 2 (l-z 2 ) 



sD, 



d=6.EWl 




(l + pz)(l-p z 



(3 - p 2 + 2pz( 1 + pz) - ( 1 + p 2 + 2pz) p z + 2p 2 ) gU 



+ 2(j3z+(l + pz)p z -p 2 j^: 
1 1 



sD\ 



d=6,EW2 




) 



0,-50* (l-Pz)C^(l,2,4)-(l + Pz)C 2 (l,2,4) +0(8), 



647T 



(11.30) 



with 
B = 



-8naa s (N 2 - 1) — ^ ((J - 2 - p 2 (l - z 2 ))^, + pz( J - 2) (J - 3)^ 

(11.31) 

Note that the contribution rfa EW ~ box is UV-finite, as already mentioned. This can be 
easily checked by replacing the Bq integrals by A and verifying that this contribution 
indeed vanishes. Since we are using box integrals in d = 6 dimensions, the IR sin- 
gularities appear only in the three-point integrals. In particular, in the above result, 
only the last line in Eq. (III. 301) is singular. As a consequence only this term needs to 
be evaluated in d dimensions. Note that owing to the structure of the IR-singularities 
in QCD we will not pick up finite terms of the form e/e. They will cancel with the 
corresponding terms in the real corrections, as we will show in the next section, where 
the real corrections are discussed. 



For the QCD-box we find 
do QCD a 1 



dz 



= c 
+ 2 
+ 



71 (l-p z ) 

1 



2 



So(l,3) 



Z ( - 1 - 2P 2 + PV) gU ~ P ( 1 + P V) gtta 



*o(M) 



l-pz 



(1 - 2P 2 + pv)^:, - pcp - z) (1 + wj a 
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+ 



1 + pz 

1-p 



(l - 2p 2 + p 2 z 2 )gU v + P(P +z) (1 - Pz)^ 



^o(2,4) 



P(l-pV 
+ p(i + p¥)^ ; 

p 2 (l-z 2 ^ 



z (l+(3 2 + (3 2 (l-(3 2 )(l-z 2 



>VOV 



QCD1 



(1,3,4) 



+ 



1+pz 
p 2 (l-z 2 ^ r 
1-Pz 



3 + 2p z (l+Pz)-p 2 W v + 2pz^ 



sD 
sD 



4=6,2 


4=6,1 




1 



1 



64 ^a,.-^p5[(l-pz)Ci(l,2,4)-(l + pz)C 2 (l,2,4)j+0(£). (11.32) 

Combining the two results in Eqs. (III.30I) and dll.321) . the IR-divergent part is thus 
given by 

OCv 1 



3271 JV 2 



Pfl (1 - Pz)Cj(l,2,4) - (1 + pz)C 2 (l,2,4) 



a, 1 P^ 



\6ilN 2 s 
In particular, we obtain 

_o^^P 
16% N 2 s 
a s 1 p 



V C^1,2,4)- V -C 2 (1,2,4) 



(11.33) 



B 



Sfl ,Cj(l,2, 4) -^(1,2,4) 
2 



327rA^ 2 s 



fl(4n)T(l+e) 



■In 



Sqt 
Sqt 



+ g(Sqt)-g(s q t) 



(11.34) 



with s 



2ki ■ kj, where we have used 



<3(1,2,4) = ^/^ 



£ 2 (£ + ^) 2 ((£ + ^) 2 -^ 2 ) 

- Wr(l +£)-3-(i + lln +g( % )) +0(8), (11.35) 



with 



We note that the divergent three-point integral 



Sqt - W 



2 



V 



1 



1 



m 2 i £ 2 (£ + ^) 2 ((£ + ^ + ^-) 2 -mz 2 )' 
which could in principle also appear, cancels in the calculation. 



(11.36) 
(11.37) 
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III. Real corrections 




Figure III. 1 : Sample diagrams for the real corrections. 



As mentioned in the previous section the contribution from the box diagrams is IR- 
divergent. To render the corrections to the total cross section finite we need to include 
the real corrections at the same order. A few sample diagrams are shown in Fig. ML 1 1 
The diagram containing the triple gluon vertex (see Fig. MI.2l) does not contribute be- 
cause of the colour structure. The calculation of the real corrections is straightforward. 
The phase-space integration over the regions where the emitted gluon is soft will pro- 
duce the IR singular contribution needed to cancel the corresponding singularities in 
the virtual corrections. Note that owing to the colour structure no collinear singularities 
appear, because the interference between the two diagrams, where the gluon is emitted 
from the initial state, vanishes. As a consequence no factorization of initial-state sin- 
gularities is required. To extract the IR divergences, we use the so-called subtraction 
method ll2Tll22ll23ll . The basic idea of the subtraction method is to add and subtract a 
term in such a way that the singularities appearing in the real corrections are matched 
point- wise and that the term is simple enough to be integrated analytically in d dimen- 
sions over the full phase space. Given that the same term is added and subtracted, this 
procedure does not change the result. The analytically integrated term is combined 
with the virtual corrections, while the unintegrated term is combined with the real cor- 
rections. Given that the term combined with the real corrections match point- wise the 
singularities of the squared matrix element, the integration can be done numerically 
in 4 dimensions. Because of the universal structure of soft and mass singularities in 
QCD, the subtraction terms can be constructed in a very general way. For further de- 
tails on the subtraction method, we refer to Refs. ET1I23II . Here we just reproduce the 
necessary equations required for the case at hand. 

Using the subtraction method the NLO contribution to the cross section can be sym- 
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Figure III.2: Amplitude containing the triple gluon vertex. The diagram does not 
contribute here because of the colour structure. 



52 {dOLO ® ^Vdipole) (k q , kq,k t ,kt, k g 

dipoles 



e=0 



bolically written as ll2"Tll2"3l 

VNL0(kq,kq) =Ov(k q ,kq) +<3 R (k q ,kq 

[do R (kq : kq,k t , kj, kg)] £=Q 

+ J [dGv(k q ,kq,k t ,ki) + da LO (k q ,kq,k t ,kt) <S> (I q + I q;]E __ {) 
+ J dx Jda L o(xkq,kq,k t ,kt)®(K q +Pq) 

+ J dxjdO L0 {kq,xkq,kt,kt)^{Kq+Pq). (III.l) 

Here Jgr, dcy denote the real and virtual corrections to the cross section. In particular 
we have 

Jo v = do lR - + do Fm - + do EW - h0X + do QCD - h0X . (EI.2) 

In Eq. Jin.lt we label the integral symbols with an index 2 or 3 to indicate that the 
phase-space integral runs over a 2- or 3-particle final state. The terms of the form 

dc LO ®F (HI.3) 

with F = P, K, JVdipole deserve some explanation. In general the symbol '®' introduces 
spin as well as colour correlation between the operator F and the leading-order am- 
plitude, which is a vector in colour space. Note that for the case studied here, where 
the gluon is always emitted from a quark line, no spin correlation appears. For the 
contribution from the integrated dipoles we obtain 

do LO (k q ,kq,k t ,kt)®l = * 4^2 ~ ~^ti\l q + lq\qq ^ tf)^ 
a, 1 |3 (Anf r2, (s 



-«(v)} 



32n 2 AN 2 s r(l-e) U 
[qq -► tt I T q T t \qq->tt), (III.4) 
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with 



His. 



2 ^ ( Sqt ) ^ ( Sqt ) ^ ( 







Q 2 (s qt ,m t ) 



- In 



m t 



Q 2 {s qt ,m t 

2 
Sqt 



In 



Sqt 



Q 2 (s qt ,m t ) 



+ ln \ s qt / \ V / 2 \"v) + 2 ln ( V 



+ In 



31n 



Q 2 (s qt ,m, 
Q(s qt ,m t ) -m t 



2Li 2 



»9* 



ln 



Q 2 (s qtl m t )J s qt \Q 2 (s qt ,mt 



3m t 



- - 

Q(s qtl m t ) +m t 6 ' 



(III.5) 



Q{s qt .,m t ) 

Following Ref. lEDl we used the bra-ket notation to represent the leading-order ampli- 
tude as vector in colour space: 

\qq — > it) . 

In the derivation we used the following result for I 



I q {e^ 2 ;{},p. 



a s (47t) £ 
27ir(l-e) 

,2 ' 



1 

C~F 



T t T q 



C F 



V 



VQ(s qt ,m t ,0;£,K) 



+ r Q (v,m t ;£)+y q ln(^-^J+y q +K q 



+ h T,T « 

+ (f «->i 



C F 



Sqt 



V q (s qt ,0,m t ;e,2/3)- 



7T 



+ ^+y g in 



Sqt 



+ Jq+K q 

(m.6) 



which can be easily obtained from Ref. [23]. The definitions of 1/q, V q ,T Q,y qi K q can 
be found in Ref. fH\. The result for 1^ can be obtained from the above result by the 
replacement (q <-> q). The T t , T q appearing in the above equation are colour-charge op- 
erators, which act on the leading-order amplitudes that are vectors in colour space. The 
calculation of this specific contribution was further simplified by noting that, because 
of the simple colour structure of the process at hand, the following relation holds: 

(qq — ► tt | T q T t | qq — > it) = — (qq — > tt | T q Tj\ qq — ► ff ) 
= — {qq-+tt\ T q T t \qq -> tt) = (qq —> tt\ T q Tj\ qq — > tt) . 

The square of the colour-correlated tree amplitudes is given by 

(qq — > tt | T q T t \qq~^tt) 

= -32n 2 aa s (N 2 - 1) — ^ ((d - 2 - (3 2 (1 - z 2 ))gH + W ~2)(d- 3)gU a 

s — ntz v 

= AnB, (III.8) 



(III.7) 
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with B as defined in Eq. (HOTTl . Comparing Eq. (HCT) with Eqs. (Himi . (UlLSll . it is 
easy to see that when combining the real and virtual corrections the IR singularities 
indeed cancel. In addition, as promised already in the previous section, we see that 
indeed no e/e terms appear, which must be due to the 'J-dimensional' factorization 
of the infrared singularities. The K and P operators can be calculated along the same 
lines as described above for the I operator. In particular we obtain 



~^Z T t T q{ K ( x ^ s qti m t) - K ( x , s qh m t ) } , 



(III.9) 



with 



K(x,Sqt,m t 



'gQ 



m t 



+ 2 



LI-jcJ 



In 



(2-x)s, 



qt 



+ 8(1-*) 



In 



m t 



1 m t 2 3 
+ ~ ; ^ + ^- 



K 2-x)s qt +m t < 
2m t 



s q t \s qt +m t 2 J 2s qt + m t 2 2 y/ / Sqt + m 2 + mt 

(1 -X)Sqt 



3 / s qt ~ 2m f a/ Sqt + m t 2 + 2m t 2 
+ -In 1 



(1 -x)s qt + m t 



where the regular part F r e| of the evolution kernel P qq (x) is given by 

Pgj s (x) = -(l+x), 
and the Plus prescription defines distributions in the usual way, through 

[F(z)] + = Km |®(1-z-ti)F(z) -8(1 -z-Ti)jT V(y)rfy}. 



(III. 10) 



(111.11) 



(III. 12) 



The function is given in Eq. (5.58) of Ref. Il2"3l . and we reproduce it here 

explicitly: 



+ 



l-JC 



+ 



2(1 -x + y 2 ) 2 l-x 
1 ln(2+/-jc) 



(l+ln(l-jc + y 2 )) 



+ 



2-x 



(III. 13) 



For the P operator we find 



\ = ^ q {x)T t T q \n 



Sqt 



with 



l-x 
14 



+ 



+ -8(l-x) 



(III. 14) 



(III. 15) 



Note that in deriving the above relations we used again relation Eq. (1111.71) to simplify 
the colour charge algebra. The corresponding results for the antiquark give the same 
contribution for the K and P operators. Note that when calculating 

J dx J^dGLo(xk q ,kq,k t ,kt)®(K q + F q ), (TH.16) 

one has to replace k q by xk q in the calculation of the colour-correlated matrix element 
Eq. (IIII.8I) as well as in the phase-space measure. For details concerning the evaluation 
of the Plus prescriptions we refer to Ref. [23 1- A remark might be in order concerning 
the appearance of the K and P operators. To the order we are working here, there is 
no contribution from the factorization of initial-state singularities. At first sight the 
fact that the evolution kernel P qq appears might thus look a bit strange. The reason 
is just that a corresponding term is included in the dipole contribution, which is com- 
bined with the real corrections. If we apply the subtraction method as it is described 
in Ref. [23 1, we thus have to consider the contributions from the K and P operators 
as shown above. In principle one could think of changing slightly the form of the 
subtraction terms; but in that case the analytic integration over the subtraction terms 
would also need to be redone. Let us close this section with some remarks about the 
subtraction term 

(JO LO ® rfVdipole) (kq, k q , k t , kj, kg) (III. 17) 

in Eq. (Mill) , which is combined with the real corrections. This contribution is ob- 
tained as a sum over individual 'dipoles' <Df-, <D l J: 



tg 



+ <D? g + <D?*+0¥ g +!Df g ). (III. 18) 



Here i and j are the unresolved partons, while a plays the role of a spectator. The ex- 
plicit expressions for the 8 dipoles can be easily obtained from Ref. Il2"3l . For example, 
we get 



,/,;<;, = \ Ly<? 

(k t +k g ) 2 -m t 2 x tg:q tg 



x \q(k q )q(k q ) t(k t )t(k t ) T q ■ T t q(k q )q(k q ) -> t(k t )t(k t ) 
1 1 



_y1 ft 

(kt + k g ) 2 - m t 2 x tm tg 



- - , (HI.19) 

kq=k q ,k t =k t 



where 



*« = ^xtti t' k ' - (IIL20) 

K = x tm k» % = % + W-(\-x tm )k>i, (111.21) 
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and 

V t q g = S%a s { -I -~z t - ^— } . (IH.22) 

The four-momentum of parton i is denoted by . The momentum fraction z f is defined 
by 

~ * 23 ) 

kq ~ kf -\- kq • kg 

For the remaining dipoles the obtained results are similar. 

The numerical implementation of the subtractions terms shown above is straightfor- 
ward. So far we have only discussed the calculation of the total cross section. In prin- 
ciple also more exclusive quantities can be calculated without any significant change. 
In the next section we will discuss our numerical results for the cross section. More 
exclusive quantities will be discussed elsewhere. 



IV. Numerical results 




Figure IV. 1 : Schematic representation of real corrections to the QCD-box contribution 
(a) and the EW-box contribution (b). 



In the following we will discuss numerical results for the weak corrections to the total 
cross section. If not stated otherwise, we used the following values for the masses 

m z = 91.1876 GeV, m w = 80.425 GeV, m H = 120 GeV, 

m b = 4.82 GeV, m t = 178.0 GeV, 
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and 

a ( 2m ') = 7^' a, = 0.1, sw 2 = 0.231 
Ho. j 

for the couplings. Using mz and as input parameters, the weak mixing angle can 
be in principle calculated within the theory. However, since our calculation is leading 
order in the electroweak coupling, we expect that the numerical choice for sw as given 
within the MS-scheme, will give results closer to the actual values. Before showing 




-3 -2-1 2 3 

10 10 10 1 10 10 10 



Figure IV.2: Comparison with results available in the literature [24J. 

our final results for the cross section, let us first discuss several checks we performed. 
The analytic expressions for the vertex corrections to the initial vertex can be found in 
the literature. We compared with the results given in Refs. [ 25 , 26 1 and found complete 
agreement. For the corrections to the final vertex, no such compact expressions can be 
found in the literature. Using the same input parameters as in Ref. lfT3l we compared 
with plots shown in Ref. [ 13 1 and found agreement. Furthermore a precise numerical 
comparison with Bernreuther, Fucker and Si [29| who also recently finished an inde- 
pendent calculation of the weak corrections, lead to complete agreement. For the case 
of the box diagrams it is possible to compare the contribution of the QCD boxes with 
an analytic result available in the literature ll2"4l . In Ref. [24] the corrections to the to- 
tal cross section were calculated using the optical theorem. In our calculation we just 
need to split the real corrections into the contribution belonging to the EW-box and that 
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-40 1 1 

-3 -2-1 2 3 

10 10 10 1 10 10 10 



Figure IV. 3: Different contributions to the electroweak corrections for incoming up- 
quarks: Initial vertices (long-dashed), final vertices (dotted), EW-box (dash-dotted), 
QCD-box (dashed). The sum is shown as a full line. 
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belonging to the QCD box. As far as the matrix elements for the real corrections are 
concerned, the contribution where the momentum flow in the Z-propagator is the total 
centre-of-mass energy belongs to the QCD box. On the other hand, the contribution 
where the momentum of the gluon propagator is equal to the total centre-of-mass en- 
ergy belongs to the EW-box contribution. Sample diagrams for both contributions are 
shown Fig. II V. 1 1 For the subtraction terms there is no such difference, they have to be 
distributed equally. This splitting might sound somewhat artificial, but it allows a di- 
rect comparison with Ref. In Fig. IIV.2l we show the analytic result from Ref. [ 24 1 
as a line. The crosses are obtained from the numerical integration of our results for the 
QCD box over the full phase space. We find complete agreement taking the numerical 
uncertainties of the phase-space integration into account. This is a highly non-trivial 
test, because the entire contribution from the subtraction method is checked. In ad- 
dition we compared again with the results by Bernreuther, Fucker and Si and found 
perfect agreement [29|. 

Note that for the box contributions only the axial-vector part (proportional to g { a gt) 
contributes to the total cross section thanks to the Furry theorem. We included the 
vector part (proportional to g' v gv) in our calculation as well, because our aim is to allow 
also the calculation of differential quantities (with or without cuts) where these terms 
might contribute. We checked that for the total cross section the vector part indeed 
cancels in the numerical evaluation of the phase-space integrals — providing a further 
check of our numerical implementation. Terms proportional to g[ga or g' a gy contribute 
to parity violating observables only. They are relevant for spin dependent quantities 
and have not been included in the present analyses. An important consequence of 
the Furry theorem is that the result for incoming down-type quarks can be obtained 
directly from the one for up-type quarks as far as the boxes are concerned. There is 
just a relative sign between the two contributions, because of the sign difference in the 
weak isospin. At the hadron level the contribution of the box diagrams is thus directly 
proportional to the difference of the parton distribution functions between up- and 
down-type quarks. This leads to a suppression of the contribution of the box diagrams. 

Let us now discuss the numerical results for the cross section. In Fig. IIV.3l we show the 
separate contributions as well as the sum for the partonic cross section for incoming 
up-quarks as a function of 

" = 5?- L <IV1) 

For a Higgs mass of m# = 120 GeV used in Fig. IIV3I the dominant contribution is 
given by the vertex corrections. It can also be seen that the contribution from EW- 
boxes is much larger than the contribution from the QCD-boxes. We have checked 
that the purely weak contributions, of order a 2 are completely negligible. At hadron 
colliders the parton subenergies may reach the order of TeV and beyond. In this re- 
gion the suppression of the cross section by large Sudakov logarithms starts to be- 
come important, similiarly to the situation for purely electroweak processes (see e.g. 
Refs. lfT5l[T6l0 . In this region the weak corrections are of the order of 10% and more. 
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Furthermore, the vertex corrections depend strongly on the Higgs mass as shown in 




Figure IV.4: Relative change of the cross section from corrections to the final vertex 
for different Higgs masses m# = 120 GeV (full line), m# = 180 GeV (dashed), m# = 
240 GeV (dotted), m H = 1000 GeV (dashed-dotted). 



Fig. irV.41 and this dependence is particularly pronounced in the threshold region. (For 
a related discussion see Ref. f27|.) In particular for small m# and small velocity [3 one 
is still sensitive to the attractive Yukawa force. For large m# and small p\ on the other 
hand, the Higgs contribution vanishes. In Figs. HV.51 and IIV61 the total contribution 
of the weak corrections to the quark-antiquark induced part of the hadronic cross sec- 
tion is shown, which should be compared to the QCD corrected total cross section of 
5.75 pb OH and 833 pb El for 1.96 TeV and 14 TeV respectivly. As far as the total 
cross section is concerned, the effects are evidently very small. Electroweak correc- 
tions are, however, important for differential distributions, which are enhanced in the 
region of large parton subenergies ll30ll . 
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Figure IV.5: Dependence of the qq induced hadronic cross section on the top mass 
m t (LHC) for 3 different Higgs masses (m# = 120 solid line, m H = 200 dashed line, 
mu — 1000 dashed-dotted line). 
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Figure IV.6: Dependence of the qq induced hadronic cross section on the top mass m t 
(Tevatron) for 3 different Higgs masses (m# = 120 solid line, m H = 200 dashed line, 
m,H — 1000 dashed-dotted line). 
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V. Conclusion 



In this article we have evaluated the complete electroweak corrections to top-quark 
pair production in quark-antiquark annihilation including terms from the interference 
between QCD and electroweak amplitudes. In particular we present short analytical 
results, which will be useful for further investigations. As a first application we have 
studied the impact of the weak corrections on the total cross section. We confirm the 
findings of Ref. 1 13 1 that the correction to the total cross section is negligible. 

Acknowledgments: We would like to thank W.Bernreuther, M. Fucker and Z.-G. Si 
for useful discussions and for a detailed comparison of results prior to publication. 



A. List of used integrals 

Using the definitions 

5 (Pl,mi 2 ,m 2 2 ) = -X ! d d £j^ 5 r^?^ 5 r 

UVFU i» in 2 J ^2_ mi 2 + i - e )(^ + /7l )2_ m2 2 + J - e ) 

Co(pi,p 2 ,Pl-P2,mi ,m 2 ,m 3 ) = 



in 2 J {I 2 - mi 2 + ie)((£ + pi) 2 - m 2 2 + ie)((£ + pi+ pi) 2 - m 3 2 + ie) 
the integrals used in section[TT]are 



4(1,3) 
4(1,4) 
4(3,4) 


= B (s,0,m z 2 ) 
= B (m 2 ,0,m t 2 ) 
= B (m t 2 ,mz 2 ,m t 2 ) 




(A.l) 
(A.2) 
(A.3) 


4(2,4) 


= 5 (-^(l-pz)+m, 2 ,0,m r 2 ) 




(A.4) 


4(2,4) 


= B (~(l + $z)+m t 2 ,0,m t 2 ) 




(A.5) 


4(1,3) 

4(1,2) 
4(1,3) 

4(1,2) 
4(1,3) 


= B (s,m t 2 ,m t 2 ) 
= BQ{m 2 ,m b 2 ,m w 2 ) 
= B (s,m b 2 1 m b 2 ) 
= Bo(m 2 ,m 2 ,m H 2 ) 
= B (s,0,0) 




(A.6) 
(A.7) 
(A.8) 
(A.9) 
(A. 10) 


0,(1,3,4) 


= C (s,m t 2 ,-^,0,m z 2 ,m t 2 ) 




(A.ll) 


Cj(l,2,4) 


= Co(0,-|(l-(3 Z )+m f 2 ,^(l-(3z; 


,0,0,m f 2 ) 


(A.12) 


Q 2 (l,2,4) 


= C (0 -*-(! + £z)+m 2 /-(l+ fc] 


,0,0,m f 2 ) 


(A.13) 
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Cj(2,3,4) = C (0,m, 2 ,-^(l-pz),0,m z 2 ,m ? 2 ) (A.14) 



Cj(2,3,4) = C (0,m f 2 ,-^(l + |3z),0,m z 2 ,m ? 2 ) (A.15) 
C ^ CD1 (1,3,4) = C (s,m f 2 ,-J,0,0,m, 2 ) (A.16) 



s 



2' 

C = Co(m t ,m t i2~ mt > m < ' mz ,m? ) ( A - 17 ) 



C = C (m f ,m, ^~ mt ltTlb ' mw ,m& - 1 (A. 18) 

C = C (m f ,m t ^2~ mt > m < > mH ,m? ) ( A -19) 

of 6,W1 = D6(0,0,m, 2 ,|,-^(l-(3z),-|(l + Pz),0,0,m z 2 ,m^.20) 

= D^(0,0,m, 2 ,|,-^(l + Pz),-|(l-Pz),0,0,m z 2 ,m^.21) 

Df 6 ' 1 = D6(0,0,m f 2 ,|,-^(l-(3z),-^(l + |3z),0,0,0,m f 2 ) (A.22) 

Z)f 6 ' 2 = D6(0,0,m f 2 ,|,-^(l + (3z),-^(l-(3z),0,0,0,m f 2 ) (A.23) 
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